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ABSTRACT 

We demonstrate a new way of studying interplanetary magnetic field - atomic alignment. In- 
stead of sending thousands of space probes, atomic alignment allows magnetic mapping with 
any ground telescope facilities equipped with spectro-polarimeter. The polarization of spectral 
lines that are pumped by the anisotropic radiation from the sun is influenced by the magnetic 
alignment, which happens for weak magnetic field (< 1G). As a result, the line polarization 
becomes an excellent tracer of the embedded magnetic field. The method is illustrated by the 
specific cases of Io and comet Halley that we consider. Magnetometer data from the Galileo 
mission Io flyby (2002) were used in order to construct the topology of the magnetic field 
around Jupiter. So as to the data from the vega mission comet Halley flyby(1986). A uniform 
density distribution of Na was considered and polarization at each point was then constructed. 
Both spatial and temporal variations of turbulent magnetic field can be traced with this tech- 
nique as well. For remote regions like the the boundary with interstellar medium, atomic 
alignment provides a unique diagnostics of magnetic field, which is crucial for understanding 
the physical processes like the IBEX ribbons discovered recently. 

Key words: interplanetary medium, magnetic fields, atomic processes, scattering, techniques: 
polarimetric 



1 INTRODUCTION 

Magnetic fields are ubiquitous and play extremely important roles in many astrophysical circumstances, e.g., the interstellar medium, inter- 
galactic medium, and quasars. In the mean time, there are only a few techniques for magnetic field studies that are available. Each technique is 
sensitive to magnetic fields in particular environments and has its own limitations. Therefore, even the directions of a magnetic field obtained 
for the same region of sky with di fferent techniques diffe r substantially. The simultaneous use of different techniques provides a possibility 
of magnetic field tomographv.tsee iYan. & Lazarianll 20061 for details) 

Interplanetary magnetic field is turbulent with complicated spacial and temporal variations. However, it is always important or unavoid- 
able when considering the interplanetary processes. For instance, we have to deal with the coherence of the plasma and the magnetic field 
when studying the interplanetary medium. The interplanetary current sheet is created by the magnetic field's rotation and reconnection. The 
comic rays can be obviously influenced by the interplanetary magnetic field. 

The recently launched NASA's satellite Interstellar Boundary Explorer (IBEX) l McComas et al. 20091) discovered a bright narrow 
"ribbon" of energetic neutral atoms (ENA) emission nearly across the entire celestial sphere, which is unpredicted by any prior models or 
theories. It has been shown that this ribbon is la rgely stable, with s mal l time variations, which implies a sustained local magnetic structure 
interacting with the original solar magnetic field i Funsten et aljyeat and Schwadro n et alj 2009h. Galactic magnetic field immediately outside 
the heliopause's special orientation might result in the "ribbon" structure. (see Heerikhuisen et al. 2010). The heliosphere is asymmetric as 
the two Voyager spacecraft detected when crossing the termination shock (see lstone et al.l2005h . indicating the significance of magnetic field 
in the dynamic process of the outer heliosphere. But as far as we know, there is no direct measurement of the magnitude and orientation of 
the magnetic field in those region. 

Atomic alignment is a new promising way of studying magnetic fields in radiation-dominated environments. In fact, this diagnostic is 
far more sensitive than those based on the Zeeman effect. The basic idea of atomic alignment is simple: consider atoms irradiated by a nearby 
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Figure 1. Left: a toy model to illustrate how atoms are aligned by anisotropic light. Atoms accumulate in the ground sublevel M = as radiation removes 
atoms from the ground states M = 1 and M = — 1 ; Right: radiation geometry and the polarization vectors in a given coordinate system. Q. is the direction of 
radiation, % is the positional angle of a linear polarization. 



star and embedded in a magnetic field. Anisotropic radiation pumps the atoms differentially from different magnetic sublevels, resulting in 
over- or underpopulations of the atomic states of various magnetic quantum numbers, M. Then the absorbed and scattered light will get 
polarized by the aligned atoms, relying on the magnetic field, especially the direction, around them. So it's an effective way to obtain the 
3D orientation of magnetic field and a cost efficient method, just by detecting the polarization of light from the object rather than sending a 
satellite. 

In this paper, we first discuss a synthetic observation of Jupiter's magnetic field that has been detected by Galileo spacecraft. The third 
largest of Jupiter's moons, Io, is chosen as the tracer of our synthetic observation, since it has been shown that there is abundant sodium in the 
wake of Io. Secondly w e give the result of synthetic observation of the Comet Halley. The alignment of sodium was calculated particularly in 



Yan. & Lazarianl 1 2007 ). Here, we use the D2 emission line of sodium {high resolution). We also consider the situation when the spectrometer 



cannot resolve the Dl and D2 lines (low resolution). Since Dl line is unpolarized, the polarization degree in the case of low spectral resolution 
is reduced. 

Galileo spacecraft was sent to Jupiter for an 8-year mission in Jovian system. It detected Io and other moons of Jupiter, se nding abundant 



inform ation about the Jovian system. One of them is the magnetic field. We downloaded the data from NASA PDS-PPI fsee lKivelson et al 



1997). 



Vega spacecrafts (vegal and vega2) were sent to Venus and to Comet Halley was implemented by the USSR with international partic- 
ipation. The main spacecraft encounter with Comet Halley was in March 6 and 9, 1986. Data were obtained during both the encounter and 
cruise phases of the missionQ 

Undoubtedly, compared to the vast interplanetary space that magnetic fields prevail, the information we can get from space crafts is 
far less adequate. With these studies, we try to demonstrate atomic alignment as an easy, but powerful magnetic diagnostics. It enables a 
tomography of interplanetary magnetic field without sending thousands of space probes. 

In what follows, we first discuss the basic idea of atomic alignment, giving some important formula that is necessary for our synthetic 
observation in §2. Then talk a little about the frame rotation in §3, which is simple but sometimes confusing. We present the observation 
result of the magnetic field traced by Io in §4 and discuss the advantage of atomic alignment to detect the interplanetary magnetic field. The 
synthetic observation of Comet Halley is discussed in §5. The discussion and the summary are provided in, respectively, §6 and §7. 



2 BRIEF INTRODUCTION OF ATOMIC ALIGNMENT 

The basic idea of the atomic alignment is quite simple. The alignment is caused by the anisotropic deposition of angular momentum from 
photons. In typical astrophysical situations the radiation flux is anisotropic. As the photon spin is along the direction of its propagation, we 
expect that atoms scattering the radiation from a light beam to be aligned. Such an alignment happens in terms of the projections of angular 



see |ftp : //nssdcftp . gsf c . nasa . gov/spacecraft_data/vega/mag/| 
2 Since the alignment of atom and ion are generally the same, atom represents both of them in this paper 
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momentum to the direction of the incoming light. For atoms to be aligned, their ground state should have nonzero angular momentum. 
Therefore, fine (or hyperfine) structure is necessary to enable various projection of atomic angular momentum to exist in their ground state. 

Let us discuss a toy model, that provides an intuitive insight in the physics of atomic alignment. Consider an atom with its ground state 
corre sponding to the total angular momentum 1=1 and the upper state corresponding to the angular momentum 1=0 (see Varshalovich et al] 



1968). If the projection of the angular momentum to the direction of the incident resonance photon beam is M, the lower state M can have 



values -1,0, and 1, while for the upper state M=0 (see FigH/e/f). The unpolarized beam contains an equal number of left and right circularly 
polarized photons whose projections on the beam direction are 1 and -1. 

Thus, absorption of the photons will induce transitions from the M=-l and 1 sublevels. However, the decay from the upper state populates 
all three sublevels of the ground state. As a result, the atoms accumulate in the M=0 ground sublevel from which no excitations are possible. 
Accordingly, the optical properties of the media (e.g. absorption and scattering) would change. 

The above toy model can also exemplify the role of collisions and magnetic field. Without collisions one may expect that all atoms 
reside eventually at the sublevel of M=0. Collisions, however, redistribute atoms t o different subl evels. Nevertheless, as disalignment of the 
ground state requires spin flips, it is less efficient than one might naively imagine i Hawkinslll955 ). The reduced sensitivity of aligned atoms 



to disorienting collisions makes the effect important for various astrophysical environments. 

All in all, in order to be aligned, first, atoms should have enough degrees of freedom; namely, the quantum angular momentum number 
must be 1. Second, the incident flux must be anisotropic. Moreover, the collisional rate should not be too high. While the latter requires 
special laboratory conditions, it is applicable to many astrophysical environments such as the interplanetary medium, the outer layers of 
stellar atmospheres, the ISM, and the intergalactic medium. 

For the optically thin case, the linear polarization degree and the positional angle are 

P=^Q 2 + U 2 /I=^ 2 2 + z 2 Jz (h x=\^-\u/Q) = ItairHea/ei) (1) 

(see Fig\T\right), where e,- (f = 1,2,3) are emission coefficients. As to the alignment of Sodium, things get somehow complicated. Since 
Sodium has nuclear spin 1 = 3/2, we sho uld consider the total mo mentum F = J + I, thus hyperfine structure is taken into consideration. We 
can obtain the emission coefficients (see lYan. & Lazarianl (2007)) for detail calculation) of each hyperfine structure as below. For D2 line, 
Fi = h 

e = !^_A/,np§£(v-v ) {62830 -10084cos2e r -412cos48 r + (30.3cos48 r +2236.6cos28 r + 776) cos 29 + cos 2<|v [(395 - 513.4cos29 r + 107.4cos49,.)sin 2 9j +cos$ r [104.3 

+ 1 .95 cos26v - 212. 19 cos48, + 8.5(cos26y - 1.4cos46, + 48.5) sin28 sin26y]} 
E] = -^^A/,np§yv-Vo) {(1552.7 + 4473.3 cos29 r + 60.6cos49 r ) sin 2 9+cos2<|>,. [592.5 -f cc*28(197.5 + 53.7 cos48 r -256.7 cos28 r ) -770 .1 cos28,+ 161.1cos48,] 

+ cos 4>,. ( 1 7 .0 cos 26,- - 23. 8 cos49 r ) sin 29 sin 29 r } 

E ; 7 = -^2^La/,kp|J4(v-Vo) [cosO (790 - 1026.9cos29 r +214.8cos49, )sin2c|>,. + (34.0cos29 r -47.6cos49 r + 1649.9) sin<|ysin9sin29 r ] (2) 

Fi = 2, 

£o = -^J— A '»"Po^( v - v o) { 136670 -28893cos29 r +899cos49,- + cos 2 8(20729cos28, - 11 18cos48 r + 6057) + cos24>,- (5649.4 - 6278.5cos29,- 

+ 656.2cos49 r ) sin 2 8 + cos<|),. [2522.6- 147.8co.s29,-- 2521.1 cos49,+ (5090 -591 cos 29,- + 6 cos 49,) sin 29 sin 29 r ]} 
£] = -^^A/»«p^(v-v ){(20729cos29,. - 1117.7cos49 r + 6057.1) sin 2 9 + cos2^) r [(-6278.5 cos29 r + 656.2cos49 r + 5649.4) cos 2 9 - 6479.6cos29,- + 656.2cos49 r + 5649.4] 



+ cos 4>,.( 10097- 1183.0cos29,.+ 12.4cos49 r )sin29sin29,. + 6057.1} 

_3\/3X. : 
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A/,np[j5(v -V ){cos9 (29189- 48 18cos29 r + 215 cos49 r ) sin 2 9,- sin2^),. + (20193 -2366cos29,. + 25cos49,.)sin<j),.sin9sin28,.| (3) 



For Dl line, F; = 1, 

E„ = l^^A/.iipgyv -v n ) {36142 - 6350cos29 r - 285cos49,- - cos29(40 + 76cos49 r + 26cos29,-) -cos24>,- [(588 - 702cos29,- + 1 26 cos 49, ) sin 2 9J - cos^,- [268 - 1.19cos29,- - 267cos49 r 

+ (537.0 -5.1 cos29 r + 2.7cos49 r ) sin29sin29,[ } 
E, = -^2^_A/,np[j5(v-V( ) ){(130.89cos29 r - 153.45cos49 r -2050.2) sin 2 9 + cos20,- [cos29(-293.92 - 62.56cos49,. + 251.03cos29,.) + 1053.08cos29 r - 187.674cos48, 

881. 772] -cos4>,. [1073.92- 10.216cos29 r + 5.472cos49 r ]sin29sin29 r + 1969.6} 
E 2 = - ^2^^A/,/ip^(v-v ) [1404. 1 cos9(cos29, -0.1782cos49 r - 0.8373) sin2<ti,. + (20.4cos29,.-10.9cos49 r -2147.8) sin^,-sin9sin29,] (4) 



En 



F,=2, 



A/»iip o ^(v-v ) {64488- 10148cos29, + 337 cos49 r -cos29(-66.73cos49, +65.45cos29,- -40.301) -cos2^),- [(-587.8 + 702cos29,- - 125.2cos49,-) sin 2 9] 



8Tt 

+ cos4>,- [268 - 1 . 19cos28 r - 267. 1 cos49 r + sm2esin2e,.(536.96 -5.11 cos29 r + 2.74cos49 r )]} 
El = -62« = 1) 

E2 = -E2(F/ = 1) (5) 
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Figure 2. Left: typical astrophysical environment where atomic alignment can happen. A pumping source deposits angular momentum to atoms in the direction 
of radiation and causes differential occupations on their ground states. In a magnetized medium where the Larmor precession rate V£ is larger than the photon 
arrival rate T^ 1 , however, atoms are realigned with respect to magnetic field. Atomic alignment is then determined by 9,-, the angle between the magnetic 
field and the pumping source. The polarization of scattered line also depends on the direction of line of sight, 9 and 9o (or § r , defined afterward). Right: The 
xyz-frame is JSO (or CSE) frame in which the data were detected. We firstly need to transmit the magnetic data and the space coordinates into observational 
frame. This can be done by two successive rotations specified by Euler angles (9,0). The first rotation is from xyz coordinate system to x'y'z' coordinate 
system by an angle <]) about the z-axis, the second is from x'y'z' coordinate system to x"y"z" coordinate system by an angle 9 about the y'-axis. 



3 FRAME ROTATION 

Since the theoretical calculation of alignment is generally in different frame from the observational data recorded in, some frame rotations 
are necessary, (see FigJT]nght, Fig^right and Fig\3\right) We get Jupiter's magnetic data from NASA PDS-PPI in two separate frames, JSO 
(Jupiter Solar Orbitaljjand SYS3 (system III) both of whose instruction can be found in the corresponding webside with the data0 The 
significant difference between the two frames is that SYS3 frame is rotating together with Jupiter, the period of which is about 9 and a half 
hours, while the JSO frame is rest in Jupiter-Sun system. So the relative position of the sun, the Jupiter and the earth can approximate to be 
quiescent, since the revolution of both the Jupiter and the earth is quite long compared with the time used for collecting the magnetic field 
data. The facility of choosing those two frames is that data in SYS3 can show us the map of Jupiter's magnetic field, even though it can 
merely give a qualitative impression. And the data in JSO frame can gives the picture of polarization corresponding to the magnetic field in 
observational frame on the earth, which we mainly employ. For the case of Comet Halley, the coordinate frame is CSE (CometoCenter Solar 
Ecliptic almost the same as JSO except that the center object is the comet. 

Some initial conditions are necessary, which are the position angle of the sun (the pumping source) and the position angle of the earth 
in JSO (or CSE) frame. Since the x-axis is defined along the radial from Jupiter(or Comet Halley) to the sun, the position angle of the sun is 
simply (0°,0°)- We choose the earth is to be at (83°, 5") in the case of Jupiter and (60°, 5°) in the case of Comet Halley. Since this is just a 
synthetic observation, the coordinates are not accurate. 

In terms of the synthetic observation, we need to obtain Q r and <]),. which are crucial for calculating polarization. Firstly we need 
to transform all magnetic field data (from JSO frame) into the observational frame (see Fig(2] right), since the calculation of emission 
coefficients are based on the frame, (see Fig\3\right) And it's worthy to mention that all the figures below that's in x-y plane corresponds 
to the observational frame here we received. Secondly, in order to obtain 9,- and (j),-, we transform the observational frame to the theoretical 
frame (Fig\T\right) which are generally various in different position subjected to the direction of magnetic field. In general, those two rotations 
are the same and can be done by two Euler rotations, as illustrated in Fig^right and Fig[3\right. Mathematically, the two rotations from 
observational frame to the theoretical frame can be fulfilled by multiplying rotation matrixes, 



cos 6 — sin 6 

1 
sin 6 cosG 



cos<])£ sin(])£ 
— sin(j)g cos<j)g 
1 



cos 6 cose 

— sin 
sin 6 cos (t 



cos6sin<|)£ —sin 6 

cos<|)£ 
sin 6 sin cosG 



(6) 



3 In JSO coordinates, X points from Jupiter to the Sun. Z is parallel to the Jupiter orbital plane upward normal. Y completes the right handed set. 

4 see |http: //pds-ppi. igpp.ucla.edu/ inventory . jsp?volume=GOMW_50 05&folder=DATA/MAG/SURVEY/REALTIME&file=ORB32_JSO_JSM| 

5 In CSE coordinates, the origin is at comet position, x-axis is from comet to the Sun, z-axis is to ecliptic northpole, (x,y)-plane is parallel to the ecliptic plane. 
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Figure 3. Left: geometry of the observational frame. In this frame, the line of sight is z axis. The incident light is in (9o,4>o) direction. Magnetic field is in 
§b) direction; Right: transformation to the "theoretical frame" where magnetic field defines z" axis. The process of rotation is the same as Fig[2] right. 
Thus, The first rotation is from xyz coordinate system to x'y'z' coordinate system by an angle §b about the z-axis, the second is from x'y'z' coordinate system 
to x"y"z" coordinate system by an angle 8g about the y'-axis. Atomic alignment and transitions are treated in the "theoretical" frame where the line of sight 
is in (9,7t) direction and the incident radiation is in (Q r ,§ r ) direction. (see Fig[T]rig/!f and Fig^/e/f) 
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Figure 4. Left: the synthetic observation result in celestial coordinate, where red thick line represents the polarization direction and blue thin line represents 
the direction of magnetic field. The unit of the coordinates is Rj=71492 km. Right: the magnetic field detected by Galileo, represented by the same data with 
the left panel but in SYS3 frame. It can give us the qualitative impression of the Jupiter's magnetic field. 



The rotations from JSO (or CSE) frame to observational frame are all the same. You can just change §s to to obtain the transform 
matrix according to Fig^right. We obtain the position angle of pumping source in theoretical frame, 

0,- = arccos(sin0cos<j)flsin9ocos(])o + sin 6 sin sin Oo sin 00 + cos0cos9) (7) 

ty r — arcsin((— sini|)B sin9ocos(])o + cos(|)b sinBo sini|)o)/sin6 r ) (8) 

where equals 9# and (0q, (|>o) is the position angle of the pumping source in observational frame, (see Fig(3j 
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Figure 5. Left: a local part of the Fig^'e/f which distinguish the blue thin lines. Right: the acute angle between the polarization vector and the projection line 
of the magnetic field vector to the celestial coordinate plane. To ensure the clarity of the figures, we use lines instead of arrows to represent the vectors. 
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Figure 6. Left: the polarization degree of low resolution. Right: the polarization degree of high resolution, n of the x-axis is the sequence number of the data 
corresponding to Fig^/e/i". The two panels correspond to FigHjLe/i. 
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Figure 7. Left: 9,-, the angle between incident radiation and magnetic field see Fig|2] /e/f and Fig[5] right. Middle: <|v, the azimuth of the incident radiation 
corresponding to 9,- in the theoretical frame. Right: 9, the angle between the observational direction and the magnetic field. 
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Figure 8. Left: the Magnetic field and polarization along the trajectory of vegal near the comet Halley. The blue thin lines represent magnetic field and the red 
thick lines represent the polarization vectors. The x-y plane is the plane of celestial coordinates. Right: the acute angle between the polarization vector and the 
projection line of the magnetic field vector to the celestial coordinate plane, n of the x-axis is the sequence number of the data, corresponding to the left panel. 



4 SYNTHETIC OBSERVATION 

Since spacecrafts can only get the local information of magnetic field on its orbit, our synthetic observation shall follow the trajectory 
of spacecrafts. Magnetic field increases sharply when it goes near Jupiter, as shown in Fig\4\right. Correspondingly, the blue thin lines are 
distinct only at the bottom of FigUfe// which locates near Jupiter, although magnetic field exist everywhere. In order to reveal the relationship 
of the magnetic field and the polarization, we chose and magnified a local part of the Fig[4]fe/f as Fig[5]fe/f. %, the acute angle between the 
polarization vectors (red thick lines) and the projection of the magnetic field in the celestial coordinate plane (blue thin lines, Fig\5\right), 
varies in the range of ±50 degree (see Yan. & Lazarian 1 20071) ). For the low resolution, the average polarization degree is around 3 %, while 



it can reach as high as 7.5 % and as low as 0.5 % (see Fig[6] left). Fig|6] right shows us the polarization degree at high resolution, whose 
polarization degrees are generally higher than low resolution s, since the light of Dl line is unpolarized. The mean polarization degree is 4.3 
%, while it can reach 11.1 % to the most and 0.8 % to the least. It's easy to find that the shapes of the polarization degree figure at high and 
low resolution is the same, the difference is the range. The similarity of Fig[7]/e/f with right is rational since there isn't much difference in 
the position angles of the sun and the earth in JSO frame. 



5 COMET HALLEY 

The mean polarization degree of low resolution is 4.5 %, while it can reach 10 % to the most and almost to the least, (see Fig(9]/e/0 The 
mean polarization degree of high resolution is 6.8 %, while it can reach 15 % to the most and almost to the least. (Fig[9]ng/i?) Since comets 
are cruising around the interplanetary space, it's efficient to detect the magnetic field of some particular area by tracing a comet. 



6 DISCUSSION 

Although we focus on the emission line of Na Dl and D2, there are abundant lines alignable in prominent universe lines. In other words, this 
method is practical to detect magnetic field of the interplanetary, interstellar and intergalactic space. It also can be used for the magnetic field 
in QSOs and other objects. As long as the magnitude of the magnetic field is sub-gauss and the collisional rate is not too high. 

Atom alignment is sensitive to the orientation of the magnetic field. We can acquire the entire information of the magnetic field with 
two or more alignable lines. 

The resolution of our synthetic observation is limited by the spacecraft's interval of detection. The space resolutions employed in the 
synthetic observation are FigH/e/f 6.3 x 10 5 km and Fig[8] left 1.01 x 10 4 km. Compared with the image a fine telescope can obtain in 
interplanetary space, th ey are too rough. For example the Steward MMT with resolution 0.018 arcsec at the wavelength of 0.5 /jm (see 
Strittmatter et al. 1 1979h ). can resolute the interval less than 100 km when take pictures around Jupiter. Hence, the real observation can be 



quite precise even more than the spacecrafts. In addition, the figure of all the parameters shown above seem much more like broken lines 
rather than curves also because the big interval between two data points. 

Since it's quite easy to translate the polarization information to the magnetic information, we can obtain the instantaneous tomography 
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Figure 10. Left: Q r , the angle between incident radiation and magnetic field see Fig|2]fe/? and Fig\3}right. Middle: <]),., the azimuth of the incident radiation 
corresponding to 9,- in the theoretical frame. Right: 9, the angle between the observational direction and the magnetic field. 

of the magnetic field, which is important to observe the evolution of particular magnetic field directly. Furthermore, it's possible to acquire 
the 3D topology of the magnetic field. 

7 SUMMERY 

1. Atoms and ions with fine or hyperfine structure^ > 1) can be aligned in their ground state, providing a promising method to detect the 
weak magnetic fields, especially the orientation, in diffused medium. 

2. The polarization degree P and the angle between the magnetic field and the polarization % are determined by the position angles of 
pumping source and the observer. 

3. The magnetic field magnitude does not apparently limit the polarization degree, insuring the method's validity for weak magnetic 

field. 

4. In general, the polarization degree P varies from to 15 %. 

5. The resolution of the spectrometer doesn't limit the accuracy of detecting the orientation of the magnetic field. 

6. The resolution of detecting magnetic field with this method is determined by the resolution of the telescope. The accuracy can be 
even higher than the space probe. 

7. Atomic alignment allow us to obtain the instantaneous tomography of the magnetic field. 
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